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We show how the transverse t + t~ spin correlations can be used to 
measure the parity of the Higgs boson and hence to distinguish a CP-even 
H boson from CP-odd A° in the future high energy accelerator experi- 
. ments. We investigate the subsequent decays of the into 7r ± z/ T (^ T ) and 

j o ± f' r (^ r ). The prospects for the measurement of the Higgs boson par- 
ity with a mass of 120 GeV are quantified for the case of e + e~ collisions 
of 500 GeV center of mass energy. The Standard Model Higgsstrahlung 
i' production process is used as an example. 
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^ ' 1. Introduction 

H ' 

In the Standard Model (SAi) of elementary particle interactions, the 
breaking of electroweak symmetry is achieved through the Higgs mecha- 
nism. The simplest realization is provided by the introduction of a complex 
Higgs doublet, which leads to the presence of a neutral CP-even Higgs boson 
in the physical spectrum. Among the possible extensions of the SAi, the 
Minimal Supersymmetric Standard Model (AiSSAi) has been considered 
most seriously. The minimal realization of the Higgs mechanism within su- 
persymmetric extensions of the standard model requires the presence of two 
Higgs doublets at low energies. After the Higgs mechanism operates, five 
real fields remain, and there should be five spin zero Higgs fields, and the 
spectrum includes also a pseudoscalar CP-odd state. This non trivial as- 
signment of the quantum numbers requires the investigation of experimental 
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opportunities to measure the parity of the Higgs boson states. The inves- 
tigation of the mechanism of the electroweak symmetry breaking is one of 
the central task of a future e + e~ linear collider operating at center-of-mass 
energies between 350 and 1000 GeV . This accelerator will allow to study 
completely and with high accuracy also the profile of the Higgs sector. 

The problem of the Higgs boson(s) parity measurement was approached 
in general form quite early in Ref. ^[2] for Higgs decays into fermions and 
gauge boson pairs. It has recently been revived [SlIlJOIlinHHEllHl . 

In this study we investigate the case of a Higgs boson which is light 
enough that the W + W~ decay channel remains closed. Then the most 
promising decay channel of neutral Higgs particles is the bb channel, B1Z{H — 
bb) ~ 90%. This applies in the SA4 as well as in the minimal super- 
symmetric extension like MSSM ^Uj. However, due to depolarization 
effects in the fragmentation process, it is very difficult to extract informa- 
tion on the b polarization state A much cleaner channel, though with 
branching ratios suppressed by an order of magnitude, is the t + t~ mode, 
B1Z{H —* t + t~) ~ 9%. The r + r _ channel is useful in the SAi for Higgs 
masses less than ~ 140 GeV. Up to this mass, the Higgs particle is very 
narrow, T(H) < 10 MeV. In Supersymmetric theories, the t + t~ channel 
is useful over a much larger mass range. The main production mechanism 
of the SM Higgs boson in e + e~ collisions in the future Linear Collider 
(LC) are the Higgsstrahlung process, e + e~ — ► ZH and the WW fusion pro- 
cess, e + e~ — > W*W* — > v e v e H. The cross section for the Higgsstrahlung 
process scales as ~ 1/s and dominates at low energies while the cross sec- 
tion for the WW fusion process rises as ~ log(s/mj I ) and dominates at 
high energies ^Uj- At y/s = 500 GeV, the Higgsstrahlung and the WW 
fusion processes have approximately the same cross section for 100 GeV 
<m H < 200 GeV. 

In our analysis, we will take as an example the e + e~ — * ZH; Z — * p + p~; 
H — ► t + t~ production process. We discuss a method for the parity mea- 
surement of the Higgs boson with a mass of 120 GeV for the case of e + e~ 
collisions with yfs = 500 GeV using H/A° — * t + t~; — > n v T (v T ) and 

— > p^v T (y T )\ — > vr^vr decay chains. All the Monte Carlo samples 
have been generated with the TAU0LA library jl2U13lll4j . For the production 
of the r lepton pairs the Monte Carlo program PYTHIA 6.1 is used [15j . 
The effects of initial state bremsstrahlung were included in the PYTHIA gen- 
eration. For the r lepton pair decay with full spin effects included in the 
H — * t + t~; — * 7r ± P T (^ T ) and r ± — > p v T {y T )\ p — > 7r ± 7r° chains, the 
interface explained in Refs. ^OE] was use d- It is an extended version of 
the standard universal interface of Refs. [181119] . 

The rest of the paper is organized as follows. In section 2 the theoretical 
considerations used to understand the decay chain of the Higgs particle are 
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explained. In section 3 and 4 we define the observables we use to distinguish 
between the scalar and pseudoscalar Higgs boson in the — > tt^v t {v t ) and 
r ± — > p^i) T (y T ) decays respectively. In these sections we list our assump- 
tions on detector effects and the imposed cuts as well as the main numerical 
results. A summary closes the paper. 

2. Higgs boson parity 

The H/A parity information must be extracted from the correlations 
between r + and r~ spin components which are further reflected in corre- 
lations between the r decay products in the plane transverse to the t + t~ 
axes. This is because the decay probability, see e.g. Ref. [S], 

T(H/A° -> r+O ~ 1 - sfsf ± s T ^sl (1) 

is sensitive to the polarization vectors s T and s T+ (defined in their 
respective rest frames, the z-axis is oriented in the r~ flight direction). 
The symbols ||/_L denote components parallel/transverse to the Higgs boson 
momentum as seen from the respective rest frames. This suggests that 
the experimentally clean t + t~ final state may be the proper instrument to 
study the parity of the Higgs boson. 

Now, a few representative examples of the r lepton decay will be dis- 
cussed in more detail. In particular, we analyze the r decay into one and 
two pions. 

3. Higgs CP from — >■ iz^Dt^Vt) decay 

Exploring transverse spin correlations in the Higgs boson decay H/A° — ► 
t + t~; — > 7r ± P r (i/ T ) one can provide a model independent parity deter- 
mination test. The method relies on the properties of the Higgs boson 
Yukawa coupling to the r lepton, which in the general case can be written 
as f(a T + i6 r 75)r. The method does not depend on the Higgs boson pro- 
duction mechanism at all. Even though the r ± — ► 7T P T (v T ) decay mode is 
rare, B1Z{t — > ttv t ) ~ 11%, it may serve as a simple example that illustrates 
the basic principles. 

The spin of the r ± leptons is not directly observable but translates 
directly into correlations among their decay products. To demonstrate this 
we define the polar angles between the 7r ± and the r~ direction in the r 
rest frames by ± and the relative azimuthal angle 0* between the decay 
planes. The angular correlation between t + t~ decay products may then be 
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written as [201 
1 dT(H/A° -» ir + v T ir-v T ) 1 



[l + cos# cos 9 + =F sin # + sin 9 cosi 



r d cos 9 + d cos 9 ~d(j)* 8tt 

A simple asymmetry in the azimuthal angle that projects out the parity 
of the particle can be derived by integrating out the polar angles 



1 dT{H/A°) _ 1 
F ddf ~~ 2vr 



7T 2 



1 =F cos ' 
lb 



(3) 



A useful observable sensitive to the parity of the decaying Higgs particle is 
the angle 5 between the two charged pions in the Higgs boson rest frame 
|2L)| . This delicate measurement will require, however, the reconstruction 
of the Higgs boson rest-frame. It is generally accepted that the Monte 
Carlo method is the only way to estimate whether the measurement can be 
realized in practice, and which features of the future detection set up may 
turn out to be crucial. To enable such studies we have extended the standard 
universal interface [TS ^ fTO j . of the TAUOLA r-lepton decay library [T2lll3U14j . 
to include the complete spin effects for r leptons originating from the spin 
zero particle |lbU17| . The interface is expected to work with any Monte 
Carlo generator providing Higgs boson production, and subsequent decay 
into a pair of r leptons. 

Let us turn to the discussion of numerical results. As an example we 
took a Higgs boson of 120 GeV. Fig. ^presents the distribution in the angle 
4>* = arccos(n + • n~) where 

\p w x p \ 

i.e. the acoplanarity angle. Thick lines will denote predictions for the scalar 
Higgs boson and thin lines for the pseudoscalar one. The distribution is 

2 

indeed, as it should be, proportional to ~ iTfg cos (/)* respectively for scalar 
and pseudoscalar Higgs boson, see Eg.(JSJ). In Fig. [5]we plot the distribution 
of the 7r + 7r~ acollinearity angle (5*). The difference between the case of 
a scalar and a pseudoscalar Higgs boson is clearly visible, especially for 
acollinearities close to ir. 

To test the feasibility of the measurement, some assumptions about the 
detector effects have to be made. We include, as the most critical for our dis- 
cussion, effects due to inaccuracies in the measurements of the tt^ momenta. 
We assume Gaussian spreads of the 'measured' quantities with respect to 
the generated ones. For charged pion momentum we assume a 0.1% spread 
on its energy and direction. 
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Fig. 1. The ir + ir~ acoplanarity distribution (angle (j>* ) in the Higgs boson rest 
frame. The thick line denotes the case of the scalar Higgs boson and thin line the 
pseudoscalar one. 
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Fig. 2. The 7r + 7r~ acollinearity distribution (angle 6* ) in the Higgs boson rest frame. 
Parts of the distribution close to the end of the spectrum; 5* ~ 7r are shown. No 
smearing is done. The thick line denotes the case of the scalar Higgs boson and the 
thin line the pseudoscalar one. 
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If the information on the beam energies and energies of all other observed 
particles (high pt initial state bremsstrahlung photons, decay products of 
Z etc.) are taken into consideration the Higgs boson rest frame can be 
reconstructed. We may define the "reconstructed" Higgs boson momentum 
as the difference of the sum of the beam energies and the momenta of all 
visible particles, that is decay products of the Z and all radiative photons 
of | cos ^| < 0.98. In our study we will mimic in a very crude way the 
beamstrahlung effects, assuming only a flat spread over the range of db 5 
GeV for the longitudinal component of the Higgs boson momentum with 
respect to the generated one 1 . 




Fig. 3. The tt + tt~ acollinearity distribution (angle S* ) in the reconstructed Higgs 
boson rest frame. All smearing is included. Parts of the distribution close to the 
end of the spectrum; 5' ~ n are shown. The thick line denotes the case of the 
scalar Higgs boson and the thin line the pseudoscalar one. 

Fig. |21 shows us the distribution of the acollinearity angle (<5*) build 
from smeared tt^ momenta defined in the reconstructed Higgs boson rest- 
frame. The difference between the scalar and the pseudoscalar Higgs bosons 
is only barely visible. We have studied several mechanisms of Higgs boson 
production, in all cases depletion of the acollinearity distribution sensitivity 
to transverse spin effect was quite similar. We can conclude that our results 
are thus independent from the production mechanism. 

The beamstrahlung effect taken into account in the reconstruction of the 
Higgs boson four-momentum degraded the method of measuring the Higgs 

The typical spread for the beam energy in a linear collider is of the order of a few 
percent |21|. 
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boson parity using the decay chain H — ► r + r~, r — > 7r v T (y T ) in a decisive 
way. Therefore, there is little hope for the elegant method to check the Higgs 
boson parity using its decay to — > t:^i> t {v t ), whatever the luminosity 
of the future Linear Collider is assumed, unless other, unfortunately less 
sensitive, to spin, than — > tt^Vt^t), decay modes are used as well. 



4. Higgs CP from — ► p^i> T (v T ) decay 

The discussion of the r decay to a two pion final state follows very 
much the same line. We may treat the hadron system as a single particle 
with definite spin and mass, m p = y/ Q 2 , or we may extract additional 
information from the individual pion momenta. The first strategy leads 
to a simple generalization of the single pion case. However, the effect is 
diminished as compared with the r ± decays to single t:^i> t {v t ). This is 
because for the — > p^Vr^Ur) decay the angular correlation term is reduced 
by the factor (m 2 - 2Q 2 ) 2 /(m 2 T + 2Q 2 ) 2 

1 dT(H/A° -> p+v T p-u T ) 



1 

8^ 



(m 2 - 2Q 2 ) 2 

1 + ~, \ ~~7rr~77 [COS COS 9 + =F SU1 9 + Slli 9 COS 

{m 2 + 2Q 2 ) 2 L 



r d cos 9 + d cos 9~d<p* 

' (5) 



The mass of the hadronic system, Q 2 , can no longer be neglected relative 
to m 2 (22]. 

Let us turn to the case when additional information coming from p^ 
decay products can be used directly [T?j- The — > p u T (u T ) decay is 
very interesting because it has, by far, the largest branching ratio, B1Z{t — > 
pu T ) ~ 25%. The polarimetric force of this channel can be improved if 
information on the p decay products, i.e. on details of the decay — > 
7r :l= 7r P r (^ r ), is used. This is of no surprise because the polarimetric vector 
is given by the formula 

K = N(2(q ■ N)q i - q 2 N l ) (6) 

where N is a normalization function, q is the difference of the tt^ and 7r° 
four-momenta and N is the four-momentum of the r neutrino (all defined 
in the r rest frame) see, e.g. \12\ . Obviously, any control on the vector q can 
be advantageous. It is of interest to note that in the r lepton rest frame, 
when m^-i = m^o is assumed, the term 

q-N=(E v ±-Ej>)m T . (7) 
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Thus, to exploit this part of the polarimetric vector, we need to have some 
handle on the difference of the ir^ and 7r° energies in their respective lep- 
tons rest frames. Otherwise, the effect of this part of the polarimetric vector 
cancels out and one is left with the part proportional to the p (equivalently 
v T ) momentum. 

Let us now discuss a new observable which we have introduced to distin- 
guish between the scalar and the pseudoscalar Higgs boson. We advocate 
the observable where we ignore the part of the polarimetric vector propor- 
tional to the p^ (equivalently v T ) momentum in the r rest frame. We rely 
only on the part of the vector due to the differences of the ir^ and 7r° mo- 
menta, which manifests the spin state of the p^. In the Higgs boson rest 
frame the p momentum represents a larger fraction of the Higgs's energy 
than the neutrino. Therefore, we abandon the reconstruction of the Higgs 
boson rest frame and instead we use the p + p~ rest frame which has the 
advantage that it is built only from directly visible decay products of the 
p + and p~ . In the rest frame of the p + p~ system we define the acoplanarity 
angle, <p*, between the two planes spanned by the immediate decay products 
(the tt^ and 7r°) of the two p's, see Fig. |1J 




Fig. 4. Definition of the p + p decay products' acoplanarity distribution angle, if* 
in the rest frame of the p + p~ pair. 

The variable (p* alone does not distinguish the scalar and psuedoscalar 
Higgs boson. To do this we must go further. The — > 7r =l= 7r P r (z/ T ) spin 
sensitivity is proportional to the energy difference of the charged and neutral 
pion (in the r rest frame) , see formula (J7J) . We have to separate events into 
two zones, C and D, 



<p' 



Fig. 5. The p + p~ decay products' acoplanarity distribution angle, ip* , in the rest 
frame of the p + p~ pair. A cut on the differences of the 7r ± 7r° energies defined in 
their respective rest frames to be of the same sign, selection y\yi > 0, is used in 
the left plot and the opposite sign, selection yvyi < 0, is used for the right plot. No 
smearing is done. Thick lines denote the case of the scalar Higgs boson and thin 
lines the pseudoscalar one. 
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Fig. 6. The p + p~ decay products' acoplanarity distribution angle, ip' , in the rest 
frame of the p + p~ pair. A cut on the differences of the 7r ± tt° energies defined 
in their respective replacement rest frames to be of the same sign, selection 
yiV2 > 0, is used in the left plot and the opposite sign, selection yiy 2 < 0, is used 
for the right plot. All smearing is included. Thick lines denote the case of the 
scalar Higgs boson and thin lines the pseudoscalar one. 
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C : 



ym > o 



D : 



ym < o 



where, 



Vl = 



2/2 = 



(8) 



E w ± and E n o are the 7r , 7r° energies in the respective r 1 * 1 rest frames. 

In Fig.|S]we plot the distribution of (p* , where the left hand plot contains 
the events where the energy difference between the ir + and tt° defined in 
the t + rest frame is of the same sign as the energy difference of tt~ and 
7T° defined in r~ rest frame (selection y\D2 > 0). The right hand plot 
contains the events with the opposite signs for the two energy differences 
(selection yiy2 < 0). It can be seen that the differences between the scalar 
and pseudoscalar Higgs boson are large. If the energy difference cut was 
not applied, we would have completely lost sensitivity to the Higgs boson 
parity. 

Unfortunately, since the r-lepton is not measurable, such a selection cut 
cannot be used directly. We reconstruct the r lepton rest frames, and we 
assume Gaussian spreads of the 'measured' quantities with respect to the 
generated ones. For charged pion momentum we assume a 0.1% spread 
on its energy and direction. For neutral pion momentum we assume an 
energy spread of ~?=|| =^ ■ F° r the 6 and (j) angular spread we assume | . 

These neutral pion resolutions can be achieved with a 15% energy error 
and a 27r/1800 direction error in the gammas resulting from the 7r° decays. 
These resolutions have been approximately verified with a SIMDET j23J, a 
parametric Monte Carlo program for a TESLA detector |24| . In Ref. ^7] 
the method of reconstruction of replacement lepton rest frames was 
proposed. We will use this method here as well. We boost the 7r + , 7r°, tt~ , ir° 
momenta to the respective replacement lepton rest frames. The ir^ 
energies defined this way are used in the y\ and yi energy difference cuts. 

When we use the selection cuts y± and yi and the replacement rest 
frames as well as smearing the ir^ momenta, we obtain the results shown in 
Fig. H3 We see that the effects to be measured diminish but remain clearly 
visible. 



We have studied the possibility of distinguishing a scalar from a pseu- 
doscalar couplings of light Higgs boson to fermions using its decay to a 
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pair of t + t~ leptons and their subsequent decays to r — ► n v T (y T ) and 
— > p £v(z/ r ); p 1 * 1 — ► vr 1 ^ . 

The beamstrahlung effect in reconstruction of the Higgs boson four- 
momentum degraded the method of measuring the Higgs boson parity using 
the decay chain H — ► t + t~, — ► n v T (v T ) in a decisive way. Therefore, 
there is little hope for the elegant method to check Higgs boson parity using 
its decay to — ► -k^v t (v t ) : whatever the luminosity of the future Linear 
Collider is assumed. 

For — > p ± u T (y T ) decay, we have discussed an observable which is 
very promising. Using reasonable assumptions about the SA4 production 
cross section and about the measurement resolutions we find that with 500 
fb^ 1 of luminosity at a 500 GeV e + e~ linear collider the CP of a 120 GeV 
Higgs boson can be measured to a confidence level greater than 95%. We 
emphasize that the technique is both model independent and independent 
of the Higgs boson production mechanism and depends only on good mea- 
surements of the Higgs boson decay products. Thus, this method may be 
applicable to other production modes including those available at proton 
colliders as well as at electron colliders. 
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